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Abstract 

In this paper, we study the effects of Generalized Uncertainty Principle(GUP) 
and Modified Dispersion Relations(MDRs) on the thermodynamics of ultra-relativistic 
particles in early universe. We show that limitations imposed by GUP and particle 
horizon on the measurement processes, lead to certain modifications of early uni- 
verse thermodynamics. 
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1 Introduction 

Generalized Uncertainty Principle is a common feature of all promising candidates of 
quantum gravity. String theory, loop quantum gravity and noncommutative geome- 
try(with deeper insight to the nature of spacetime at Planck scale), all indicate the modifi- 
cation of standard Heisenberg principle [1-10]. Recently it has been indicated that within 
quantum gravity scenarios, a modification of dispersion relation(relation between energy 
and momentum of a given particle) is unavoidable[ll-13]. There are some conceptual rela- 
tions between GUP and MDRs. These possible relations have been studied recently [14, 15]. 
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These quantum gravity effects, in spite of being small, are important since they can modify 
experimental results. There are several efforts to provide experimental evidence of these 
small effects. For example, Amelino-Camelia et al, by investigation of potential sensitivity 
of Gamma-Ray Burster observations to wave dispersion in vacuo, have outlined aspects 
of an observational programme that could address possible detection of these quantum 
gravity effects[16]. Amelino-Camelia and Piran have argued that Planck-scale deformation 
of Lorentz symmetry can be a solution to the Ultra High Energy Cosmic Rays(UHECR) 
with energies above the GZK threshold and the TeV-7 paradoxes [17]. Gambini and Pullin 
have studied light propagation in the picture of semi-classical spacetime that emerges in 
canonical quantum gravity in the loop representation [18]. They have argued that in such 
a picture, where space-time exhibits a polymer-like structure at microscales, it is natural 
to expect departures from the perfect non-dispersiveness of ordinary vacuum. They have 
evaluated these departures by computing the modifications to Maxwell's equations due 
to quantum gravity, and showing that under certain circumstances, non-vanishing correc- 
tions appear that depend on the helicity of propagating waves. These effects could lead to 
observable cosmological predictions of the discrete nature of quantum spacetime. Then, 
they have addressed to observations of non-dispersiveness in the spectra of gamma-ray 
bursts at various energies to constrain the type of semi-classical state that describes the 
universe. Jacobson et al have shown that threshold effects and Planck scale Lorentz viola- 
tion are combined constraints from high energy astrophysics [19]. These literatures provide 
possible experimental schemes for detection of small quantum gravity effects. However, 
there are two extreme domains: black hole structure and early stages of the universe 
evolution where these quantum gravity effects are dominant. Corrections to black hole 
thermodynamics due to quantum gravitational effects of minimal length and GUP have 
been studied extensively(see [20] and references therein). On the other hand, part of the 
thermodynamical implications of GUP and MDR have been studied by Amelino-Camelia 
et al[21] and Nozari et aZ[22]. Thermodynamics of early universe within standard Heisen- 
berg principle has been studied by Rahvar et al[23]. Since quantum gravitational effects 
are very important in early stages of the universe evolution, it is natural to investigate 
early universe thermodynamics within GUP and MDRs frameworks. Here we are going 
to formulate thermodynamics of ultra-relativistic particles in early universe within GUP 
and MDRs frameworks. In the first step, using GUP as our primary input, we calculate 
thermodynamical properties of ultra-relativistic particles in early universe. In formula- 
tion of the early universe thermodynamics within GUP framework, due to limitations 
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imposed on the measurement processes, two main points should be considered: first due 
to casual structure of spacetime, maximum distance for causal relation is particle horizon 
radius and secondly, there is a minimum momentum imposed by GUP which restricts 
the minimum value of energy In the next step, for a general gaseous system composed 
of ultra-relativistic particles, we find density of states using MDRs with Bose-Einstein or 
Fermi-Dirac statistics and then thermodynamics of the system will be followed. In each 
step we discuss ordinary limits of our equations and we compare consequences of two 
approaches. 

2 Preliminaries 

Emergence of the generalized uncertainty principle can be motivated and finds support 
in the direct analysis of any quantum gravity scenario. This means that GUP itself is a 
model independent concept. Generally, GUP can be written as [24] 

5x5p> -(l + k(5x) 2 + r]{5p) 2 + 7), (1) 

where k, rj and 7 are positive and independent of 5x and Sp (but may in general depend 
on the expectation values of x and p). This GUP leads to a nonzero minimal uncertainty 
in both position and momentum for positive k and r)[24\. If we set k — we find 

5x5p>^(l + V (5p) 2 + 1 ). (2) 

Since we are going to deal with absolutely smallest uncertainties, we set 7 = from now 
on. So we find 

5x8p>^(l + r,(5p) 2 ). (3) 

This relation leads to a nonzero minimal observable length of the order of Planck length, 
(Sx) min = hyfij. Any position measurement in quantum gravity has at least (Sx) min as 
its lower limit of position uncertainty. This relation has an immediate consequence for 
the rest of statistical mechanics: it modifies the fundamental volume u of accessible 
phase space for representative points. In ordinary statistical mechanics, it is impossible 
to define the position of a representative point in the phase space of the given system 
more accurately than the situation which is given by {Sq5p) m i n > %. In another words, 
around any point (q,p) of the (two dimensional) phase space, there exists an area of the 
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order h which the position of the representative point cannot be pin-pointed. In ordinary 
statistical mechanics we have the following definition of fundamental volume 



u = (5q5p) m . (4) 
Since in quantum gravity era 5p ~ p, we can interpret equation (3) as a generalization of 

h, 

h eff = h(i + VP 2 ). 



Therefore, we find the following generalization of the fundamental volume 

K) e// = [Mi+w 2 )r = (^//) 3jv . (5) 

Since the total number of microstates is given by Q = — (here u is the volume of the 
accessible phase space), we see that GUP leads to a reduction of accessible microstates 
and therefore a reduction of entropy. In other words, when we approaches Planck scale 
regime with high energy and momentum particles, the volume of the fundamental cell 
increases in such away that eventually the number of microstates tends to unity and 
therefore entropy vanishes. This is a novel prediction of quantum gravity. Recently we 
have calculated microcanonical entropy of an ideal gaseous system and we have observed 
an unusual thermodynamics of systems in very short distances or equivalently very high 
energy regime [22]. 

Another consequence of GUP in the form of relation (3), has been formulated by Kempf 
et a/[24]. They have shown that within the momentum representation, the generalization 
of the scalar products reads 

J 1 + T]p z 

where and ip are momentum space state functions. For ultra-relativistic particles with 
E = pc, we should consider the following generalization 

where we have set c = 1. 

On the other hand, if we set i] = in (1), we find 

5x5p>^(l + K(5x) 2 + (8) 
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where for positive k leads to nonzero minimal uncertainty in momentum. This statement 
leads to a space-dependent generalization of H. This type of generalization has nothing to 
do with dynamics and there is no explicit physical interpretation of it at least up to now. 
From another perspective, in scenarios which consider spacetime foam intuition in the 
study of quantum gravity phenomena, emergence of modified dispersion relations takes 
place naturally [25]. As a consequence, wave dispersion in the spacetime foam might 
resemble wave dispersion in other media. Since Planck length fundamentally set the 
minimum allowed value for wavelengths, a modified dispersion relation can also be fa- 
vored. Recently it has been shown that a modified energy-momentum dispersion relation 
can also be introduced as an observer-independent law[26]. In this case, the Planckian 
minimum-wavelength hypothesis can be introduced as a physical law valid in every frame. 
Therefore, the analysis of some quantum-gravity scenarios has shown some explicit mech- 
anisms for the emergence of modified dispersion relations. For example, in the framework 
of noncommutative geometry and loop quantum gravity approaches this modified disper- 
sion relations have been motivated(see for example[21] and references therein). In most 
cases one is led to consider a dispersion relation of the type (note that from now on we set 
c = h — k,B — 1) 

(p) 2 = f(E, m- l p ) ~ E 2 - /i 2 + ai l p E 3 + a 2 l 2 p E A + O (l 3 p E 5 ) (9) 

where / is the function that gives the exact dispersion relation, and on the right-hand 
side we have assumed the applicability of a Taylor-series expansion for E <C l/l p . The 
coefficients ctj can take different values in different quantum-gravity proposals. Note that 
m is the rest energy of the particle and the mass parameter \x on the right-hand side is 
directly related to the rest energy, but ^ m if the oti do not all vanish. Since we are 
working in Planck regime where the rest mass is much smaller than the particle kinetic 
energy, there is no risk of confusing between m and [i. While in the parametrization of (3) 
we have included a possible correction term suppressed only by one power of the Planck 
length, in GUP such a linear-in-/ p term is assumed not to be present. For the MDR a large 
number of alternative formulations, including some with the linear-in-/ p term, are being 
considered, as they find support in different approaches to the quantum-gravity problem, 
whereas all the discussions of a GUP assume that the leading-order correction should 
be proportional to the square of l p (as has been indicated by Amelino-Camelia et al[21], 
linear-in-/p term in MDR has no support in string theory analysis of black holes entropy- 
area relation and therefore it seems that this term should not be present in MDR. Recently 
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we have shown that coefficients of all odd power of E in MDR should be zero [15]). 
Within quantum field theory, the relation between particle localization and its energy 
is given by E > ^, where 5x is particle position uncertainty. It is obvious that due 
to both GUP and MDR this relation should be modified. In a simple analysis based 
on the familiar derivation of the relation E > ^[27], one can obtain the corresponding 
generalized relation. Since we need this generalization in forthcoming arguments, we give 
a brief outline of its derivation here. We focus on the case of a particle of mass M at 
rest, whose position is being measured by a procedure involving a collision with a photon 
of energy E 1 and momentum p 7 . According to Heisenberg's uncertainty principle, in 
order to measure the particle position with precision 8x, one should use a photon with 
momentum uncertainty 5p^ > 4-. Following the standard argument [28] , one takes this 
8p-y > 4- relation and converts it into the relation 8E 1 > using the special relativistic 
dispersion relation. Finally 5E-y > 4- is converted into the relation M > -4- because the 
measurement procedure requires SE < M, in order to ensure that the relevant energy 
uncertainties are not large enough to allow the production of additional copies of the 
particle whose position is being measured. If indeed our quantum-gravity scenario hosts 
a Planck-scale modification of the dispersion relation of the form (9) then clearly the 
relation between £p 7 and SE^ should be re-written as follows 



This relation will modify density of states for statistical systems. Note that one can use 
GUP to find such relation between 5p 7 and 5i? 7 [15]. 



Now we are going to calculate thermodynamical properties of ultra-relativistic particles 
in early universe, using the generalized uncertainty principle. We consider the following 
GUP as our primary input, 



where £ is a dimensionless constant. Consider the early stages of the universe evolution. 
Analogue to a particle inside a box, in the case of the early universe one can consider a 
causal box (i.e. particle horizon) which any observer in the universe has to do measure- 
ments within this scale[29]. In the language of wave mechanics, if \1/ denotes the wave 




(10) 



3 GUP and Early Universe Thermodynamics 




(ii) 
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function of a given particle, the probability of finding this particle by an observer outside 
its horizon is zero, i. e. \^(x > horizon)] 2 = 0. From the theory of relativity, measure- 
ment of a stick length can be done by sending simultaneous signals to the observer from 
the two endpoints, where for the scales larger than the causal size, those signals need more 
than the age of the universe to be received. Looking back to the history of the universe, 
the particle horizon after the Planck era grows as if -1 , but inflates to a huge size by the 
beginning of inflationary epoch. Here H is Hubble parameter. In the pre-inflationary 
epoch, the maximum uncertainty in the location of a particle, 8x = H^ 1 results in an 
uncertainty in the momentum of the particle which is given by 

Sp>^ + nH. (12) 
This leads to a minimum uncertainty in momentum as 

(Sp)min = 7777 + nH. (13) 
Therefore, we can conclude that (assuming that p ~ 5p) 

Pmin = J^+TTH, (14) 

which leads to 

E mm = V3(^ + nH), (15) 
for ultra-relativistic particles in three space dimensions. Now, suppose that 

E n = n&, (16) 

where •d is given by 

*=J + *ff. (17) 

To obtain complete thermodynamics of the system, we calculate partition function of 
the system and then we use standard thermodynamical relations. In classical statistical 
mechanics, partition function for a system composed of ultra-relativistic noninteracting 
monatomic particles (Fermions or Bosons) is given by 

roc AttTI 2 

lnZ = ±gj o -^ln(l±e- /3En )dn. (18) 
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In our case, due to limitation imposed by GUP and particle horizon, we should consider 
the following generalization 



InZ = ± 



gn 



E 2 



: 2^k„iT^ ul ^ e ' > ur " (19) 

where we have used relations (7), (15) and (16) respectively. By definition, the entropy 
of the system is given by 



where F is the free energy of the system defined as 



F — In Z. 



(21) 



So the entropy of the system can be written as 

E 2 . . ,„ gn(3 r°° 

' If 1 



*4 



± 



gn r°° 
2ffi J E min 1 



E 3 dE 
+ V E 2 e? E ± 1 



(22) 



For ultra-relativistic fermions this relation leads to the following expression 



S f = 



7 7T 4 31 777T 6 1016?7 2 7r 8 3 (l r]E, 

— — - -— 7Z 1 — — E m -\ J 



7 2tt 2 (1 + D) [90 /3 3 210 /3 5 1680 /5 



77 2 £ 4 . 
mm i / mm 



3' ° 5 

where for simplicity we have defined 



-^/ 3 + ^/3/ 5 -^ 2 /3/ 7 + ... 



7 J ln(l+e _/3E ™") 
(23) 



■E min E^dE 



While for bosons we find 
9 



S h = 



4 7T 4 48 r]n 6 64 r] 2 n 8 



2vr 2 (l + D) [45/5 3 315 /? 5 105 /3 7 



min 



\ V^min 1 ^ ^nm j ln(l — e"^" 1 ™) 



-^J 3 + ^J 5 -^ 2 /?J 7 + ... 



(24) 



where 



J 4 



rE m i n E^dE 
ePE-l' 
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In these equations D is denned as 

D 



,A 3 A 2 A, 



and d = A + B 



with A = jfjj and B = ttH. Note that both of the equations (23) and (24) are well 
behavior in high and low temperature limits. In the standard situation, we have £ = 0, 
n = and E m i n = 0. So we find the well-known and standard results for entropy of 
the corresponding ultra-relativistic fermionic or bosonic systems. From (22) we find the 
following expression for standard entropy 



S = 



which leads to 



and 



4 (3g r°° E 3 dE 
32W0 e^±l' 

9 7tt 4 



2tt 2 90/? 3 ' 
9 4 vr 4 



(25) 

(26) 
(27) 



2tt 2 45/5 3 ' 
for fermions and bosons respectively. 

Now the pressure of the ultra-relativistic gas is given by P = jylnZ. For fermions and 
bosons we find respectively 



Pr = 



2tt 2 (1 + D)(3 



7 vr 4 31 nn 6 127 n 2 7i 8 , . 
3^^~T260^ + T680~^ 



3 - I 

mm \ 



5 



E 2 E 4 " 

mm 1 2 mm 



+7] 



P J \ P J P 2r 



and 
P b = 



2tt 2 (1 + D)(3 



1 ^ 8 W \ 8 ^tt 8 3 (1 E 2 mm 



E 4 " 

2 min 



45 /? 3 315 /? 5 105 [3 7 



In the standard situation, we find the following well-known result 

g r°° E 3 dE 



6tt 2 Jo eP E ±V 



ln(l+e _/3£mm ) 
(28) 

(29) 
(30) 
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which leads to 



and 



2tt 2 360 p 4 ' 

p =jlL*A 

b 2tt 2 45 / 3 4 ' 

for fermions and bosons respectively. 

The specific heat of the system which is defined as 

(dS\ 



(31) 
(32) 



C v = T[ — \ 



(33) 



v 



can be written in the following closed form 

g(3 2 r™ E 4 



c 



V 



2tt 2 (1 + D) 



roc 



dE 



+ V E 2 e-P E {eP E ±lY 



(34) 



One can obtain explicit form of Cy for fermions and bosons using relation (33), (23) and 
(24). A simple calculation gives 



Cy f = 



21 7r 4 155??7r 6 7112?? 2 7r 8 
2tt 2 (1 + D) [90/33 ~~ 210^ + 1680 ~f r ~ e? E ™™ 



-4 

min 



2p4 * 

min 



+ 1\3 

8 2 dl 7 

7^ ~dT "' 



(35) 



and 

Cyb 
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2ir 2 (l + D) 



12 vr 4 240 777T 6 448^ ^ m 



7^ 



2 

min 



45 f3 3 315 /? 5 105 /3 7 e^™» -1\ 3 



+ 



r? 2 E 4 . * 

/ mm 



6 



4dJ, 6 dJ R 8 rfJ 7 



for specific heat of fermions and bosons respectively. In the standard case we find 

Cyf = 3 x — — = 3Sf, 



and 



Vb 



3 x 



2tt 2 90 f3 3 

9 4 7T 4 



3S h 



(36) 



(37) 



(38) 



2tt 2 45 [3 3 

Figure 1 shows the values of entropy in different situations. In standard thermodynamics 
of ultra-relativistic fermionic or bosonic gas, the entropy of the system tends to zero in 
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T = 0. This situation is shown in Figure 1, (a) and (b). Within GUP framework, entropy 
tends to zero in a nonzero temperature, that is, for T > T . This is a result of quantum 
fluctuation of spacetime itself. Figure 2 shows the corresponding behavior of pressure as 
a function of temperature. Note that these figures are plotted in arbitrary units and they 
show only general behaviors of the functions. Figure 3 shows the behavior of specific heat 
of the system in various conditions. In GUP framework, the general behavior of Cy has 
considerable departure from its standard counterpart in high temperature regime. 



4 MDR and Early Universe Thermodynamics 



Now we are going to formulate early universe thermodynamics within MDR framework. 
We consider a gaseous system composed of ultra-relativistic monatomic, non-interacting 
particles. First we derive the density of states. Consider a cubical box with edges of length 
L (and volume V = L 3 ) consisting black body radiation(photons). The wavelengths of 
the photons are subject to the boundary condition j- — jjr, where n is a positive integer. 
This condition implies, assuming that the de Broglie relation is left unchanged, that the 
photons have (space-)momenta that take values p — Thus momentum space is divided 
into cells of volume V p = (yr) = ft?. From this point, it follows that the number of modes 



with momentum in the interval [p, p + dp] is given by 

g(p)dp = 8irVp 2 dp 



(39) 



Assuming a MDR of the type parameterized in (9) one then finds that (m — for photons) 



and 



dp- 1 <M.h> E ~ - ^\)l 2 p E 2 



dE 



Using this relation in (39), one obtains 
g(E)dE = 8nV 



1 + 1a x l v E + 5(^a 2 + \(4)l 2 p E< 



E 2 dE. 



(40) 
(41) 

(42) 



This is density of states which we use in our calculations. Note that we have not set 
«i = to ensure generality of our discussions, but we will discuss corresponding situation 
at the end of our calculations. 
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To obtain thermodynamics of the system under consideration, we start with the partition 
function of fermions and bosons, 

In Z = ± / g(E) ln(l ± e~ pE )dE, (43) 

where + and — stand for fermions and bosons respectively and (3 = 1/T since ks — 1- 
Using equation(42) in the following form 

g(E)dE = 8vrU(l + aE + bE 2 )E 2 dE, (44) 



where for simplicity we have defined a = 2a±l p and b = 5(^a 2 + ^a 2 jll, one can compute 
the integral of equation(43) to find the following expression for entropy of fermions and 
bosons 

This relation can be written as follows 

4 j;- 9(E)H1± ^ )dE (40) 

By calculating this integral, we find for fermions and bosons respectively 



8tt 



7 tt 4 225aC(5) 31 6vr 6 „ /l a£ mi „ bE 2 -\, . RF . 

I A E • — I — A lnfl + e - ^ 

90/3 3+ 8 /3 4 + 210/3 5 rom l3 + 4 + 5 J { ' 



4 5 6 

--f31 3 --(3al 4 --(3bh + ... 
o 4 o 



and 

5 6 = 8tt 



4 tt 4 30aC(5) 48 far 6 3 /l aE miB y ^ 

45 /3 3 + /3 4 + 315 (3 5 + mm l 3 + 4 + 5 ) { 



4 5 6 

-~PJs- - A (3aJ 4 - -pbJ 5 + .. 
3 4 5 



(47) 



(48) 



One may ask about the relation between these two results and corresponding results of 
GUP, that is, relations (23) and (24). Although these results seem to be different in 
their /3 dependence, but note that if we set ol\ = (which is reasonable regarding the 
argument presented in page 5), we find a = and then (3 dependence of our findings will 
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coincide with each other. The only difference which remains is the differences between 
numerical factors. This argument shows that essentially the results of GUP and MDRs for 
thermodynamics of the early universe do not differ with each other in their temperature 
dependence and overall behaviors. 

In the standard situation, we have a = b = and E min = 0, so we find 



4 r c 



E 3 dE 
e? E ±l 

For entropy of fermions and bosons we find respectively 

7 7T 4 



and 



8tt 



8tt 



45 p 3 ' 



In the presence of MDR, the pressure of corresponding systems are 



8tt 
J 



7 7T 4 45aC(5) 31 fevr 6 



360 /5 3 8 (3 4 1260 (3 5 



1 , E n 
3 + °-4 



(49) 



(50) 



(51) 



ln(l + e"^""™) 



■3/3-^/4-56/5 + ... 



(52) 



8tt 
J 



1 tt 4 aC(5) 8 &7T 6 /I 
45 /? 3 /3 4 315 /? 5 ™"l3 4 



5 



^ 7 Put 
o 4 o 



(53) 



for fermions and bosons respectively. In the standard situation we find the following 
well-known relation 

3 Jo eP E ±r 1 ' 

which for fermions and bosons leads to 



8tt 



and 



8tt 



360 (3 4 ' 

1 7T 4 



45 P 4 



(55) 



(56) 
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respectively. 

The specific heat of the system can be written in the following closed form 



P 2 



C v = — 



g(E)E 2 dE 



V J Emin e-(> E (e?> E ± I f 



(57) 



One can use relations (33), (47) and (48) to find the following explicit results for fermions 
and bosons respectively 



C 



Vf 



8vr 



21 7T 4 900aC(5) 155 &?r e 



3E: 



rain 



90 p 3 8 p A 210 p 5 eP E ™ + l\Z 



1 aE v 



+ 



bEl 



and 



C 



Vb 



JrT ^R T ^RhT Adh 5 dh V^a 

+ 3^3 + -Jah + -fibh -a— - -b— + ... 



8tt 



PR 



•4 

rain 



12 vr 4 120aC(5) 240fe7r b 
45^ + W 4 + 315 p 5 ^ e^™" - 1 I 3 



„E ■ hE 2 ■ 



(58) 



4 5„ T 6„, T 4dJ 3 5 dJ 4 6,rfJ 5 



In the standard case we find 



90 p 3 



C vf = 3 x 8tt— — = 35/, 



and 



V6 



4 vr 4 

3 x ftr- ? = 3S b 



(59) 



(60) 



(61) 



respectively. 

As has been indicated, there are severe constraints on the functional form of MDR which 
these constraint are motivated when one compares black hole entropy-area relation in 
different points of view[15,21]. In this case we should set cti — which leads to a = 0. 
We find from (47) and (48) the following expressions for entropy of fermions and bosons 
respectively 



8vr 



90/5 3 210 p 5 mm \ 3 5 



ln(l + e" 



-0E„ 



4 6 
--/3h - -=0h + 

O 



(62) 
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and 



8vr 




45 /3 3 + 315 f3 5 1 ' ; 

4 6 
--(3J 3 --(3b'J 5 + ... 
o o 



(63) 



where b' = ^a 2 lp- These statements for partition function are more realistic since black 
hole thermodynamics within MDRs when is compared with exact solution of string the- 
ory, suggest the vanishing of a±. 

It is important to note that the formalism presented in this section is not restricted to 
early universe. Actually, it can be applied to any statistical system composed of ultra- 
relativistic monatomic noninteracting particles which has a minimum accessible energy . 
The Possible relation between GUP and MDRs itself is under investigation[14,15]. Gener- 
ally these two features of quantum gravity scenarios are not equivalent, but as Hossenfelder 
has shown, they can be related to each other[14](see also [15]). As a result, it is natural to 
expect that under special circumstances, our results for early universe thermodynamics 
within GUP and MDRs should transform to each other. This is a transformation be- 
tween coefficients of our equations and overall behaviors of thermodynamical quantities, 
specially their temperature dependence are similar. 



5 Summary 

GUP and MDRs have found strong supports from string theory, noncommutative geom- 
etry and loop quantum gravity. There are many implications, originated from GUP and 
MDRs, for the rest of the physics. From a statistical mechanics point of view, GUP 
changes the volume of the fundamental cell of the phase space in a momentum dependent 
manner. On the other hand, MDR leads to a modification of density of states. These 
quantum gravity features have novel implications for statistical properties of thermody- 
namical systems. Here we have studied thermodynamics of early universe within both 
GUP and MDRs. We have considered early universe as a statistical system composed of 
ultra-relativistic particles. Since both particle horizon distance and GUP impose severe 
constraint on measurement processes, the statistical mechanics of the system should be 
modified to contain these constraint. Since GUP and MDRs are quantum gravitation 
effects, the modified thermodynamics within GUP and MDRs tends to standard ther- 
modynamics in classical limits. There are severe constraints on the functional form of 
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MDRs from string theory considerations. When we consider these constraints, the results 
of MDRs and GUP for thermodynamics of early universe tends to each other in their 
general temperature dependence and they differs only in their numerical factors. This 
fact may be interpreted so that GUP and MDRs essentially are not different concepts of 
quantum gravity proposal. Although the exact relation between GUP and MDRs is not 
known yet, our formalism of early universe shows the very close relation between these 
two aspects of quantum gravity. 

In standard statistical mechanics of bosonic and fermionic gases, the entropy of the sys- 
tem tends to zero in To = 0. As our equations and corresponding numerical result show, 
within GUP framework entropy of the system tends to zero in a temperature larger than 
zero( T > T ). This is a consequence of the relation (5). The volume of the fundamen- 
tal cell of phase space increases due to GUP. Note that MDRs give entropy-temperature 
relation which has no difference with GUP result in its general behavior. Figure 2 shows 
the pressure of the system versus temperature. Pressure tends to zero in a tempera- 
ture larger than T = 0. The same behavior is repeated by specific heat of the system. 
So, our analysis shows an unusual thermodynamics for statistical systems in quantum 
gravity eras. This unusual behaviors have been seen in other context such as black hole 
thermodynamics[30,31] . 
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Figure 1: Entropy of ultra-relativistic monoatomic gaseous system for (a) standard bosonic gas 
(b) standard fermionic gas (c) bosonic gas within GUP and (d) fermionic gas within GUP. 
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Figure 2: Pressure of ultra-relativistic monoatomic gaseous system for (a)standard bosonic gas 
(b) standard fermionic gas (c) bosonic and fermionic gas within GUP. The difference between 
bosonic and fermionic gasses in this case is not considerable. 
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Figure 3: Heat Capacity of ultra-relativistic monoatomic gaseous system for (a)standard bosonic 
gas (b) standard fermionic gas (c) bosonic gas within GUP and (d) fermionic gas within GUP. 
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